Independent Measurement of the Total Active 8 B Solar Neutrino Flux Using an Array of 3 He 
Proportional Counters at the Sudbury Neutrino Observatory 
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The Sudbury Neutrino Observatory (SNO) used an array of 3 He proportional counters to measure the rate 
of neutral-current interactions in heavy water and precisely determined the total active (v. v ) 8 B solar neutrino 
flux. This technique is independent of previous methods employed by SNO. The total flux is found to be 
5.54^3| (stat)^Q34 (syst) x 10 6 cnT 2 s -1 , in agreement with previous measurements and standard solar models. 
A global analysis of solar and reactor neutrino results yields Am 2 = 7.59*^^ x 10~ 5 eV 2 and 6 = 34.4+] ?, 
degrees. The uncertainty on the mixing angle has been reduced from SNO's previous results. 

PACS numbers: 26.65. +t, 13.15. +g, 14.60.Pq, 95.85.Ry 



The Sudbury Neutrino Observatory Jit] detects 8 B solar 
neutrinos through three reactions: charged-current interac- 
tions (CC) on deuterons, in which only electron neutrinos 



participate; neutrino-electron elastic scattering (ES), which 
are dominated by contributions from electron neutrinos; and 
neutral-current (NC) disintegration of the deuteron by neutri- 
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nos, which has equal sensitivity to all active neutrino flavors. 

In its first phase of operation, SNO measured the NC rate 
by observing neutron captures on deuterons and found that 
a Standard-Electro-Weak-Model description with an undis- 
torted 8 B neutrino spectrum and CC, NC, and ES rates due 
solely to v e interactions was rejected at 5.3cr (2, IH 0, Etl • The 
second phase of SNO measured the rates and spectra after the 
addition of -2000 kg of NaCl to the 10 6 kg of heavy water 
(D2O). This enhanced the neutron detection efficiency and the 
ability to statistically separate the NC and CC signals, and re- 
sulted in significant improvement in the accuracy of the mea- 
sured v e and v x fluxes without any assumption about the en- 
ergy dependence of the neutrino flavor transformation 00]. 
In the present measurement, the NC signal neutrons were pre- 
dominantly detected by an array of 3 He proportional coun- 
ters (Neutral Current Detection, or NCD, array (8J) consist- 
ing of 36 "strings" of counters that were deployed in the 
D2O. Four additional strings filled with 4 He were insensitive 
to the neutron signals and were used to study backgrounds. 
Cherenkov light signals from CC, NC, and ES reactions were 
still recorded by the photomultiplier tube (PMT) array, though 
the rate of such NC events from 2 H(n,y) 3 H reactions was sig- 
nificantly suppressed due to neutron absorption in the 3 He 
strings. As described in this Letter, the new measurements 
of the CC, NC, and ES rates result in reduced correlation be- 
tween the fluxes and improvement in the mixing angle uncer- 
tainty. 

The data presented here were recorded between November 
27, 2004 and November 28, 2006, totaling 385.17 live days. 
The number of raw NCD triggers was 1,417,81 1 and the data 
set was reduced to 91,636 NCD events after data reduction de- 
scribed in ||S[] . Six strings filled with 3 He were excluded from 
the analysis due to various defects. The number of raw PMT 
triggers was 146,431,347 with 2381 PMT events passing data 
reduction and analysis selection requirements similar to those 
in Jj§|. Background events arising from (3-y decays were re- 
duced by selecting events with reconstructed electron effective 
kinetic energies > 6.0 MeV and reconstructed vertices within 
fliit < 550 cm. 

Thermal neutron capture on the 3 He in the proportional 
counters results in the creation of a proton-triton pair with a 
total kinetic energy of 764 keV. Because of particles hitting 
the counter walls [8], the detected ionization energy was be- 
tween 191 and 764 keV. The signals from each string were 
amplified logarithmically to provide sufficient dynamic range 
before they were digitized [8]. The detectors were constructed 
from ultra-pure nickel produced by a chemical vapor deposi- 
tion process to minimize internal radioactivity. 

The neutron detection efficiency and response of the PMT 
and NCD arrays have been determined with a variety of neu- 
tron calibration sources. Neutron point sources ( 252 Cf and 
241 AmBe) were frequently deployed throughout the detector 
volume to measure the temporal stability and the detector gain 
of the NCD array. The NC neutron detection efficiency was 
studied by using an isotropic source of neutrons produced 
by mixing 24 Na (ti/2 = 14.959 hours), in the form of acti- 
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FIG. 1: NCD energy spectrum fit with a neutron calibration spec- 
trum, neutron backgrounds, alpha background derived from Monte 
Carlo simulation, and low-energy instrumental background distribu- 
tions. Data are shown after data reduction up to 1.4 MeV, and the fit 
is above 0.4 MeV. 



vated NaCl, into the heavy water in October 2005 and October 
2006. Neutrons were produced by deuterium photodisintegra- 
tion induced by the 2.754-MeV 24 Na gammas. The largest 
uncertainties on the neutron detection efficiency were asso- 
ciated with the knowledge of the 24 Na source strength and 
the ability to determine the uniformity of its mixing in the 
heavy water. The inferred NC neutron capture efficiency for 
the NCD array was 0.21 1+0.007 in good agreement with the 
0.210 + 0.003 given by a Monte Carlo simulation verified 
against point-source data. The fraction of detected neutrons 
inside the analysis energy range from 0.4-1 .4 MeV, including 
the effects of data reduction, electronic thresholds and effi- 
ciency, and digitizer live time, was 0.862±0.004. The neutron 
detection efficiency for the PMT array was 0.0485±0.0006 de- 
termined from neutron point sources. 

The energy spectrum of the reduced NCD data set is shown 
in Fig.Q] The distinctive neutron spectrum peaks at 764 keV. 
This spectrum was fit with a neutron energy spectrum taken 
from the 24 Na calibration. The alpha background distribu- 
tion was derived from a Monte Carlo simulation of the pro- 
portional counters. The alpha background energy spectrum 
has several components, U and Th progeny in the bulk of 
the nickel detector bodies and 210 Po on the inner surfaces JH]. 
These sources resulted in approximately 16 alphas per day de- 
tected in the full neutron energy window for the entire NCD 
array. The Monte Carlo simulation was verified using al- 
pha data from the array above 1.2 MeV and from the 4 He 
strings in the neutron energy region. Several uncertainties 
were included in the alpha background distribution: depth 
profile and composition of alpha emitters in nickel, electron 
drift time, space-charge model parameters, and ion mobility. 
Low-energy instrumental background events were found on 
two strings that were excluded from the analysis. Distribu- 
tions of these events were used to fit for possible additional 
contamination in the data on the rest of the array. 

The optical and energy responses and position and direc- 



3 



tional reconstruction of the PMT array were updated to in- 
clude light shadowing and reflection from the NCD array. 
With improvements to the calibration data analysis and in- 
creased high voltage on the PMT array, the introduction of 
the NCD array did not significantly increase the position or 
energy reconstruction uncertainties from previous phases. A 
normalization for the photon detection efficiency based on 16 N 
calibration data [9] and Monte Carlo simulations was used to 
set the absolute energy scale. The energy response for elec- 
trons can be characterized by a Gaussian function with resolu- 
tion o- T = -0.2955 + 0.5031 y/T e + 0.0228 T e , where T e is the 
electron kinetic energy in MeV. The energy scale uncertainty 
was 1.1%. 

Backgrounds are summarized in Table Q] Low levels of 
214 Bi and 208 T1 present in the heavy and light water, NCD 
counters, and cables can create free neutrons from deuteron 
photodisintegration and low-energy Cherenkov events from 
/3 - y decays. Techniques to determine these backgrounds in 
the water are described in previous works OH U , 12 ] . In 
addition, alphas from Ra progeny on the NCD tube surfaces 
can induce 17 ' s O(a;,n) interactions. The background contribu- 
tions from the NCD array were determined by combining the 
analyses of the alpha energy spectrum and the time-correlated 
alpha events in the decay chains. The results from these stud- 
ies agreed with those from radioassays of the materials prior 
to the construction of the NCD array. In addition, in situ anal- 
ysis of the Cherenkov light found three detectable "hotspots" 
of elevated radioactivity on two strings. Evaluations of their 
isotopic composition were made upon removal of the NCD 
array after the end of data taking. Results from the in situ 
and ex situ methods showed the neutron background uncer- 
tainty from the hotspots to be less than 0.7% of the NC signal. 
Neutron backgrounds from atmospheric neutrino interactions 
and 238 U fission were estimated with NUANCE 111311 and from 
event multiplicities. 

Previous results J^, 0] reported the presence of external- 
source neutrons from the acrylic vessel and light water. Alpha 
radioactivity measurements of the acrylic vessel neck with Si 
counters before and after the NCD phase showed values con- 
sistent with those from the previous phase within uncertain- 
ties. Thus, the external-source neutron contribution from the 
vessel was taken to be the same as for the previous phase. The 
contributions from the light water were determined from the 
measured 2I4 Bi and 208 T1 concentrations. 

Backgrounds from Cherenkov events inside and outside the 
fiducial volume were estimated using calibration source data, 
measured activities, Monte Carlo calculations, and controlled 
injections of Rn into the detector. These backgrounds were 
found to be small above the analysis energy threshold and 
within the fiducial volume, and were included as an additional 
uncertainty on the flux measurements. Previous phases identi- 
fied isotropic acrylic vessel background (IAVB) events, which 
can be limited to 0.3 remaining IAVB events (68% CL) after 
data reduction for this phase. 

A blind analysis procedure was used to minimize the pos- 
sibility of introducing biases. The data set used during the 



TABLE I: Background events for the PMT and NCD arrays, respec- 
tively. Backgrounds with similar detection efficiencies are listed to- 
gether. The internal and external neutrons and the y-ray backgrounds 
are constrained in the analysis. "Other backgrounds" include terres- 
trial vs, reactor vs, spontaneous fission, cosmogenics, CNO vs, and 
(a,n) reactions. The last two entries are included in the systematic 
uncertainty estimates for the PMT array. 



Source 


PMT Events 


NCD 


Events 


D 2 radioactivity 


7.6 ± 1.2 


28.7 


+ 4.7 


NCD bulk/ 17 0, ls O 


4 6 +21 


27 6 +12 - 9 


Atmospheric v/' 6 N 


24.7 ± 4.6 


13.6 


+ 2.7 


"Other backgrounds" 


0.7 ±0.1 


2.3 ± 


0.3 


NCD hotspots 


17.7 ± 1.8 


64.4 


± 6.4 


NCD cables 


1.1 ± 1.0 


8.0 ± 


5.2 


Total internal neutron background 


56 4 +i " 


144 6 +u - 8 


External-source neutrons 


20.6 ± 10.4 


40.9 


±20.6 


Cherenkov events from fi-y decays 








IAVB 


< 0.3 (68% CL) 







development of the analysis procedures excluded a hidden 
fraction of the final data set and included an admixture of 
neutron events from muon interactions. The blindness con- 
straints were removed after all analysis procedures, param- 
eters, and backgrounds were finalized. A simultaneous fit 
was made for the number of NC events detected by the 
NCDs, the numbers of NC, CC and ES events detected by the 
PMTs, as well as the numbers of background events of various 
types. A Markov-chain Monte Carlo (MCMC) method with 
the Metropolis-Hastings algorithm lfl4l [Til was employed to 
make the fit, which also allowed nuisance parameters (system- 
atics) weighted by external constraints to vary in the fit. The 
NCD event energy spectrum was fit with an alpha background 
distribution, a neutron calibration spectrum, expected neutron 
backgrounds, and two instrumental background event distri- 
butions. The PMT events were fit in reconstructed energy, the 
cosine of the event direction relative to the vector from the sun 
(cos 8 Q ), and the reconstructed radial position. 

The spectral distributions of the ES and CC events were not 
constrained to the 8 B shape, but were extracted from the data. 



Fits to the data yielded the following number of events: 983*1? 
NC (NCD array), 267+ 24 NC (PMT array), 1867;% CC, and 
171^2 ES, with 185^2 an d 77 -m neutron background events 
in the NCD and PMT arrays, respectively. Additionally, the 
total NCD array background fits including alphas and the two 
instrumental components yielded 6127 + 101 events. 

Assuming the 8 B neutrino spectrum from lfl6ll . the equiv- 
alent neutrino fluxes derived from the fitted CC, ES, and NC 



events are (in units of 10 6 cm 2 s ') IU8U19I1 : 



,SNO 

v'cc 

,SNO 
0ES 

,SNO 
0NC 



1 67+ 005 

1 ° ' -0.04 



(stat)!°;°?(syst) 



= 1.77!0 : 2 l(stat)!» : ( »(syst) 
= 5.54-;; :; (stat)!^ (syst) 



and the ratio of the B neutrino flux measured with the CC 
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TABLE II: Sources of systematic uncertainties on NC, CC, and ES 
fluxes. The total error differs from the individual errors added in 
quadrature due to correlations. 



Source 



NC uncert. CC uncert. ES uncert. 





(%) 


(%) 


\ /u j 




+0.6 


+2.7 


+3.6 


PMT energy resolution 


±0.1 


±0.1 


±0.3 


PMT radial scaling 


±0.1 


±2.7 


±2.7 


PMT angular resolution 


±0.0 


±0.2 


±2.2 


PMT radial energy dep. 


±0.0 


±0.9 


±0.9 


Background neutrons 


±2.3 


±0.6 


±0.7 


Neutron capture 


±3.3 


±0.4 


±0.5 


Cherenkov/AV backgrounds 


±0.0 


±0.3 


±0.3 


NCD instrumentals 


±1.6 


±0.2 


±0.2 


NCD energy scale 


±0.5 


±0.1 


±0.1 


NCD energy resolution 


±2.7 


±0.3 


±0.3 


NCD alpha systematics 


±2.7 


±0.3 


±0.4 


PMT data cleaning 


±0.0 


±0.3 


±0.3 


Total experimental uncertainty 


±6.5 


±4.0 


±4.9 


Cross section f 171 


±1.1 


±1.2 


±0.5 



and NC reactions is 



jSNO 
y NC 



= 0.301 ± 0.033 (total). 



The contributions to the systematic uncertainties on the de- 
rived fluxes are shown in Table [TT] 

Two independent analysis methods were used as checks of 
the MCMC method. Both used maximum likelihood fits but 
handled the systematics differently. A comparison of results 
from these three analysis methods after the blindness condi- 
tions had been removed revealed two issues. A 10% differ- 
ence between the NC flux uncertainties was found, and subse- 
quent investigation revealed incorrect input parameters in two 
methods. After the inputs were corrected, the errors agreed 
and there was no change in the fitted central values. However, 
the ES flux from the MCMC fit was 0.5cr lower than from 
the other two analyses. This difference was found to be from 
the use of an inappropriate algorithm to fit the peak of the ES 
posterior distributions. After a better algorithm was imple- 
mented, the ES flux agreed with the results from the other two 
analyses. 

The ES flux presented here is 2.2cr lower than that found by 
Super-Kamiokande-I Ir20ll consistent with a downward statis- 
tical fluctuation in the ES signal, as evidenced in the shortfall 
of signals near cos 6 B = 1 in two isolated energy bins. The 
8 B spectral shape 111611 used here differs from that 112 ill used 
in previous SNO results. The CC, ES and NC flux results in 
this Letter are in agreement (p = 32.8% B22I0 with the NC flux 
result of the D2O phase J^t] and with the fluxes from the salt 
phase 021]. 

The fluxes presented here, combined with day and night en- 
ergy spectra from the pure D2O and salt phases |4, 7], place 
constraints on neutrino flavor mixing parameters. Two-flavor 
active neutrino oscillation models are used to predict the CC, 
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FIG. 2: Neutrino-oscillation contours, (a) SNO only: D 2 & salt 
day and night spectra, NCD phase fluxes. The best-fit point is Am 2 = 
4.57 x IO" 5 eV 2 , tan 2 9 = 0.447, f„ = 0.900, with ^ 2 /d.o.f.=73.77/72. 
(b) Solar Global: SNO, SK, CI, Ga, Borexino. The best-fit point is 
Am 2 = 4.90 x 10~ 5 eV 2 , tan 2 8 = 0.437, f B = 0.916. (c) Solar 
Global + KamLAND. The best-fit point is Am 2 = 7.59 x 10~ 5 eV 2 , 
tan 2 = 0.468, f B = 0.864. 



NC, and ES rates in SNO Q. A combined^ 2 fit to SNO 
D2O, salt, and NCD-phase data lr24ll yields the allowed re- 
gions in Am 2 and tan 2 8 shown in Fig. 12a). In a global 
analysis of all solar neutrino data (including Borexino lr25ll 
and Super-Kamiokande-I ll20Tl ) and the 2881 ton-year Kam- 
LAND reactor anti-neutrino results lE^ . the allowed regions 
are shown in Fig. [2jb and c). The best-fit point to the solar 
global plus KamLAND data yields Am 2 = 7.59!^xl0" 5 eV 2 
and 0= 34.4+JI de grees, where the errors reflect marginalized 
1-0" ranges. In our analyses, the ratio fs of the total 8 B flux to 
the SSM Ir27ll value was a free parameter, while the total hep 
flux was fixed at 7.93 x 10 3 cirr 2 s" 1 lEil . 

In summary, we have precisely measured the total flux of 
active 8 B neutrinos from the sun independently from our pre- 
vious methods. The flux is in agreement with standard solar 
model calculations. This Letter presents analysis leading to a 
reduction in the uncertainty of 9 over our previous results. 
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